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SUMMARY 

Some  approaches  using  least  squares  optimization  for  improving  tracking  control 
system  performance  are  developed  and  applied  to  a  simplified  aircraft  terrain  following 
problem. 

A  basic  tracking  control  system  is  described.  The  design  of  an  optimal  moving- 
average  precompensator  {7]  which  gives  improved  tracking  performance  over  the  basic 
system  is  then  given.  This  filter  design  is  modified  to  use  future  values  of  the  reference 
input  to  give  further  improved  performance.  An  approach  involving  the  design  of  an 
optimal  signal  to  drive  the  control  system  is  also  given. 

These  approaches  are  applied  to  an  aircraft  terrain  following  system  simulation. 
The  performance  obtained  is  examined  and  discussed.  It  appears  that  worthwhile 
performance  improvements  can  be  obtained  by  using  the  algorithms  which  make  use  of 
knowledge  of  future  terrain  as  could  be  obtained  from  a  terrain  data-base.  These 
improvements  allow  the  aircraft  to  fly  at  lower  altitudes. 

Some  proposed  extensions  of  the  work  are  described. 
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INTRODUCTION 


Of  interest  to  Flight  Management  Group  in  Aircraft  Systems  Division  at  the  Aero¬ 
nautical  Research  Laboratory  are  methods  for  incorporating  terrain  databases  into 
aircraft  involved  in  terrain  following  roles. 

In  terrain  following,  it  is  desirable  [1|  that  the  designed  aircraft  flight  path  be 
as  close  as  practical  to  the  ground  and  be  such  that  the  aircraft  is  able  to  follow  it. 
In  engineering  terms,  there  is  a  trade-off  between  following  a  desired  trajectory  and 
performing  excessive  manoeuvring,  which  includes  excessive  control  actuator  activity 
and  aircraft  normal  acceleration. 

Often,  terrain  following  system  design  treats  the  generation  and  tracking  of  the 
trajectory  as  separate  problems.  One  approach  is  to  attempt  to  follow  a  precomputed 
flight  path  consisting  of  smoothly  connected  cubic  splines,  c.g.  [1,2].  The  design  of 
such  a  path  incorporates  curvature  constraints  which  correspond  with  allowable  aircraft 
normal  accelerations.  The  path  design  may  be  formulated  as  an  optimization  problem, 
minimizing  ground  clearance  in  some  sense.  A  disadvantage  of  the  approach  is  that  not 
much  information  about  the  dynamics  of  the  controlled  aircraft  is  incorporated  into  the 
path  design. 

Another  approach  is  to  use  dynamic  programming  to  select  an  optimal  path  from 
a  prespecified  grid  of  possible  paths.  An  example  of  dynamic  programming  applied  to 
terrain  following/terrain  avoidance  is  [3].  Depending  on  the  optimality  criterion  used, 
information  about  the  aircraft  dynamics  may  be  included  in  the  path  design.  Actual 
aircraft  control  signals  may  be  generated  using  this  approach. 

A  third  approach  is  to  treat  the  terrain  following  problem  as  a  control  system 
tracking  problem  and  to  draw  from  the  large  body  of  knowledge  of  control  system 
theory  to  design  the  system.  The  problem  of  trading  desired  system  response  against 
control  activity  has  long  been  considered.  The  formulation  used  may  be  to  track  a 
precomputed  trajectory  as  in  [2]  where  a  linear  quadratic  regulator  (LQR)  is  used  to 


(1) 


track  a  cubic  .saline  path  In  (4l,  a  predictive  controller  is  applied  to  terrain  following. 
With  these  and  many  other  possible  formulations,  the  system  closed  loop  poles  result 
from  the  terrain  following  control  law  design.  There  is  then  the  possibility  that  due  to 
approximations  in  modelling  the  aircraft  dynamics,  the  system  could  be  unstable. 

The  approach  used  in  this  report  is  to  take  a  linear  dynamic  model  of  an  aircraft 
which  already  has  a  stabilizing  controller,  and  to  use  linear  control  theory  to  improve  the 
tracking  performance  while  following  a  terrain  profile,  without  affecting  the  closed  loop 
poles.  This  is  achieved  by  using  feedforward  control.  Optimal  feedforward  controllers 
have  been  described  by  Maybeck  [5]  and  Halyo  [6).  Their  work  assumes  that  the 
reference  input  is  obtained  from  the  output  of  a  linear  dynamic  system  (reference 
generator)  driven  by  white  noise.  By  solving  an  LQR  problem,  optimal  feedforward 
gains  which  operate  on  the  reference  generator  states  are  obtained. 

In  (7),  a  moving-average  (MA)  precompensator  whose  coefficients  are  chosen  to 
minimize  an  infinite  horizon  quadratic  tracking  cost  function  was  described.  The  pre- 
compensator  operates  on  actual  reference  input  values  (rather  than  on  states  of  a 
reference  generator). 

In  this  report,  that  work  is  applied  to  the  tracking  of  a  sequence  of  terrain  data 
points. 

It  is  known  [8]  that  prior  knowledge  of  the  required  plant  output  trajectory  may 
be  used  to  reduce  the  effects  of  plant  transport  delay  and  to  reduce  actuator  activity 
by  allowing  the  controller  to  respond  before  the  plant  output  is  required  to  change. 
The  terrain  following  problem  has  this  feature  if  cither  a  terrain  data  base  or  forward 
looking  radar  or  infra  red  system  is  used.  The  precompensator  design  may  be  modified 
to  make  use  of  future  reference  inputs. 

The  layout  of  this  report  is  as  follows.  In  Section  2,  the  tracking  problem  is 
formulated.  In  Section  3,  the  MA  precompensator  using  current  values  of  the  reference 
signal  [7]  is  described.  It  is  then  modified  to  allow  a  window  of  future  desired  output 


(2) 


values  to  he  operated  on  by  the  filter.  In  Section  4,  an  optimal  feedforward  control  signal 
which  is  not  constrained  to  be  a  MA  filtered  version  of  the  desired  path  is  described. 
In  Section  5,  a  simulation  example  using  an  F-15  fighter  model  and  a  section  of  real 
terrain  data  is  presented  using  each  approach. 

2  PROBLEM  FORMULATION 

In  this  work,  a  polynomial  of  order  n,  is  of  the  form 

i=0 

where  z_1  is  the  unit  delay  operator  or  the  z-transform  complex  variable.  The  argument 
(z-1)  may  be  dropped  for  brevity.  Signals  are  denoted  by  lower  case  letters,  for  example 
v(t)  or  v  for  brevity.  All  signals  considered  here  are  sampled,  with  t  =  ...,—1,0,1,... 
with  v(t)  =  0  for  i  <  0.  The  z-transform  of  v(t )  is  denoted  V(z~1 ).  The  inner  product 
of  two  square-summable  sampled  signals  r(f)  and  w(i)  with  z-tramforms  V(z  1 )  and 
\V(z~1)  is  denoted  (V,  W)  and  is  given  by 

CO 

(V,W)  =  ] Pw . 

1=0 

Matrices  are  upper  case  bold,  for  example  X,  while  vectors  are  lower  case  bold, 
for  example  x. 

The  single  input  u(t)  single  output  y(i)  plant  considered  is  described  by 

A(z~')y(l)  =  B(z~')u(t)  (1) 

where  o0  =  1  and  60, . . .  ,bq  =  0  (q  <  nj). 

This  is  controlled  by  a  general  linear  output  feedback  controller  given  by 

F(z-')u(t)  =  v(t)-G(z-')y(t)  (2) 

with  fo  =  1  where  v(t)  is  the  control  input  signal.  O  and  F  are  solutions  of  the 
Diophantine  equation 

A(z-')F(z~l)  +  B{z  ')G(z  1 )  =  T( :  ')•  (S) 


(3) 


The  closed  loop  poles  are  given  by  the  roots  of  T(z  ')  which  for  this  work  must  be 
strictly  stable  ( T”(  r  1 )  =  0  =s  |r|  <  1). 

A  schematic  diagram  of  the  system  is  shown  in  Figure  1. 

G(z~')  and  F(z~i)  do  not  appear  further  and  stabilization  via  state  feedback 
could  be  asssumed. 

Either  way, 

«(0  =  T(Z-~1\V^  ^ 

and 

m  =  Y^Mi)-  (5) 

The  desired  plant  output  signal  (also  called  the  reference  signal)  is  y"(i).  It  is 
assumed  that  ym(t)  =  0  for  t  >  ny.. 

The  cost  function  to  be  minimized  is 

oo  oo  oo 

j  =  X>*(0  -  v(<)]2  +  +  Eip^''  M*)!2-  t6> 

1=0  1=0  1=0 

The  first  term  of  J  penalizes  deviations  of  the  plant  output  y  from  its  desired  value 
y’ .  A  low  value  of  this  term  corresponds  with  accurate  tracking. 

Usually,  accuracy  of  tracking  must  be  traded  against  excessive  control  activity.  The 
second  term  of  J  in  eqn(6)  penalizes  such  excessive  activity.  Pu(z~J)  is  a  polynomial 
which  the  designer  selects  in  order  to  shape  the  frequency  spectrum  of  u(t).  If  Pu(z~l ) 
has  high  gain  at  high  frequency,  then  this  will  penalize  rapid  variations  in  u(t). 

The  third  term  of  J  allows  frequency  shaping  of  y(t). 

In  the  work  here,  methods  of  designing  v(t)  to  minimize  J  of  eqn(6)  will  be  exam¬ 
ined. 

3  MOVING  AVERAGE  PRECOMPENSATOR 

A  simple  approach  often  used  in  tracking  controllers  is  to  put 

!’(<)  -  Ky‘(i)  (7) 


(4) 


wheie  A'  is  a  constant  chosen  to  make  the  dc  gain  of  Y/Y ’  unity.  This  gives  asymptotic 
tracking  of  constant  reference  signals  and  may  also  provide  satisfactory  performance 
for  slowly  changing  reference  signals. 

The  approach  is  optimal  (minimizes  J  of  eqn  (6))  in  that  any  other  choice  of  A 
will,  for  a  constant,  non-zero  reference,  give  an  infinite  J . 

The  only  information  needed  to  calculate  A'  is  the  dc  gain  of  Y/Y’.  No  use  is 
made  of  knowledge  of  the  dynamics  of  the  system  or  the  nature  of  the  reference  signal. 

The  design  of  optimal  MA  filters  which  generate  t>  from  y*  makes  use  of  this 
additional  information. 

3.1  Moving-Average  Precompensator  Using  Current  Data 

The  controller  structure  here  is  given  by  eqn(2)  with 

«*0  =  S(z-')y’{t)  (8) 


and  the  problem  posed  is  to  design  the  precompensator  S(z~' )  to  minimize  J  of  eqn(6). 
The  designer  is  required  to  specify  weighting  polynomials  )  and  Pyi:~' )  as  well 

as  ns,  the  order  of  S(z-1). 

With  the  type  of  problem  considered  in  this  report,  for  a  given  P„  and  Py,  increas¬ 
ing  n,  gives  decreasing  J . 

The  representation  of  the  reference  as  a  polynomial  in  this  report  allows  the  de¬ 
velopment  given  here  to  be  a  simplified,  less  general  one  than  in  [7]  which  considers  the 
tracking  of  rational,  possibly  unstable  references.  Using  eHns(li)  and  (8). 


F*  -  Y  = 


Y’T  -  Y’SB 

T 


(9) 


Substituting  eqns(9),  (4)  and  (5)  into  eqn(6)  and  then  equating  the  partial  derivatives 
of  J  with  respect  to  the  s,  (i  =  0, . .  . ,  n,)  to  zero,  one  obtains  the  following  set  of  linear 
simultaneous  equations: 

Xs  -  in  (10) 


(5) 


where  s  -  (.to  •»!  ■  •  •  *«. )‘  ■ 

X  -  -  0,1,. 


VH  V'*W:  >'  Jl  ,P«Y'A  PuY'A;-''-A  /’  V*/?  IKY'liz  "  J 

A„  =■  -jO-l  +  (— f-  '  -  f - ■)  '  i  -  ,■  •  -  , 

and  m  =  (rn0  m,  ...  iiin<)J 
with  tji,  =  (^:',r*). 

X  has  a  special  structure  called  an  autocorrelation  structure  and  has  only  (u,  +  1 ) 
distinct  elements.  The  Levinson  recursion  is  an  algorithm  for  solving  equations  with 
this  structure  and  is  described  in  [9].  The  design  of  computationally  efficient  methods 
for  the  evaluation  of  the  inner  products  is  a  topic  still  of  interest  to  researchers,  for 
example  (10),  but  such  methods  were  not  used  for  this  work. 

3.2  Moving-Average  Precompensator  Using  Future  Data 

A  feature  of  the  terrain  following  problem  is  that  future  values  of  terrain  height  may 
be  available.  The  precompensator  design  is  readily  modified  so  that  the  precompensator 
operates  on  a  window  which  includes  future  values. 

Use  of  future  values  up  to  k  steps  ahead  is  achieved  by  setting 

u(f)  =  S(I-,)y*(/  +  fc)  (111 

and  the  problem  posed  is  to  design  the  precompensator  5(r~ 1 )  to  minimize  J  of  eqn(fi). 
The  designer  chooses  n„  k,  P„  and  Py.  Using  eqns  (5)  and  (11), 

.  Y‘T  -  zkY'SB 


Multiplying  by  z  k  one  obtains 


k(Y*  ~Y)  = 


hY’T  -  Y'SB 


From  the  assumption  of  zero  initial  conditions  on  y *  and  v, 


£>•(<  -  *)  -  »(*  -  *)!’  =  &•(/)  -  .v(0]2- 


(6) 


The  expression  oil  the  right  hand  side  of  eqn(!3)  can  then  be  substituted,  along  with 
equations  (4)  and  (!>),  into  eqn  (G).  Proceeding  as  in  Section  3.1.  one  obtains  the  same 
solution  (10)  but  with 

»>.  iy'y  *>'•). 

Naturally,  the  solutions  are  identical  if  k  -  0. 

4  PREPROGRAMMED  TRACKING 


Since  all  future  terrain  values  are  assumed  known,  there  is  no  need  to  constrain 
p(f)  to  be  a  MA  filtered  version  of  ym(t)  as  in  Section  3.  In  this  Section,  the  approach 
taken  is  to  design  an  unconstrained  optimal 

The  controller  structure  is  again  given  by  eqn(2)  and  the  problem  posed  is  to 
design  the  control  input  signal  v(t)  to  minimize  J  of  eqn(G).  The  designer  chooses 
n,,,  Pu  and  Py.  In  order  to  maintain  tracking  over  the  duration  of  n,.  must  be 

approximately  equal  to  ny . .  Eqn(5)  gives 


This  is  then  substituted,  together  with  eqns  (4)  and  (5)  into  eqn((i).  Partial  differ- 
entiation  of  eqn(6)  with  respect  to  the  v,  (i  -  0,...,n,,)  and  equation  of  the  partial 
derivatives  to  zero  gives 

Xv  -  m  (In) 

where  v  —  (t’o  tq  ...  i’„t  )7  , 

X  =  0,1,. 

with 


X,j  -  ( 


H  Iiz  h  ->1 
i"  T 


i  + 


Pu  A  PUA:  11  P 

r  ’  r 


Pyll  PyHz  I'  " 

T  '  T 


anil  m  ( TI10  ill 1  ...  ni„r  ) 1 
with  m,  (  H:f.  , 


(7) 


Again,  X  has  an  autocorrelation  structure.  The  Levinson  recursion  was  used  to 
solve  eqn(15)  for  v  since  nr  may  he  quite  large  («.  '  n !/-  ).  Values  greater  than  600 
have  been  used  in  this  work. 

5  SIMULATION  EXAMPLE 

The  aircraft  model  used  is  a  linearized,  discretized  model  of  the  longitudinal  air 
frame  dynamics  of  an  F-15  fighter  with  an  LQR  controller  and  was  derived  from  the 
model  used  by  Murphy  [11]  and  bv  Hill  |12j. 

The  terrain  data  used  were  a  set  of  600  collinear  terrain  height  values  obtained 
from  Woodend  in  Victoria.  These  were  spaced  93  metres  apart  ;  this  spacing  subtends 
an  angle  of  three  arc  seconds  at  the  Earth’s  centre. 

A  constant  horizontal  velocity  of  0.8  Mach  was  assumed. 

The  sampling  interval  used  for  the  model  and  the  control  update  was  0.368  seconds, 
which  is  the  time  taken  to  travel  93  metres  at  0.8  Mach.  A  new  terrain  height  value 
was  then  available  at  every  control  sample. 

The  assumption  of  constant  horizontal  velocity  causes  range  to  be  proportional  to 
time,  so  that  plots  of  performance  against  time  and  range  are  similar. 

The  plant  input,  u(f),  is  the  elevator  deflection  in  degrees  from  trim. 

Plant  output,  y(t),  is  the  aircraft  altitude  in  metres  relative  to  trim  altitude. 

The  transfer  functions  relating  u  to  v  and  y  to  e  are  given  by  eqns(  l)  and  (5)  with 

A(z“ 1 )  =  1  -  3.2452c'1  +  4.0035 c'2  -  2.5298c'3  +  1.0296c'4  -  0.2581c'5, 

=  -  0.051717c'1  -  0.738983c'2  t  0.180935c'3  t  0.564720c'4  t-  0.025510c  \ 

and 

7’(  c  '  1 )  -  1  2.4348c  1  +  2.0926c  2  0  78008c  3  t  0.16851c  4  0.043739c  5. 

The  aircraft  model  is  non-minimum  phase  with  a  zero  (root  of  li\:  1 ) )  at  c 


(8) 


8.896. 


The  desired  trajectory  y'(t)  was  obtained  from  the  terrain  profile  as  follows: 

1.  The  first  value  of  terrain  height  was  subtracted  from  all  values  so  that,  zero 
initial  conditions  corresponded  with  level  flight  at  the  initial  terrain  height. 

2.  Extra  values  decaying  the  600’t.h  element  to  zero  with  factor  0.8  were  appended 
to  }r*(c_1).  Without  this,  the  infinite  horizon  design  assumes  a  step  change  to  a  zero 
value  at  t  —  000  The  smooth  transition  to  zero  is  to  ensure  that  this  end  does  not 
affect  the  solution.  To  include  these  extra  values  ny-  =  630  was  used. 

In  order  to  penalize  excessive  differential  elevator  activity  Pu  was  chosen  to  be 

10(1 -z-1). 

Doubly  differentiated  height  is  used  to  approximate  normal  acceleration  for  small 
deviations  from  horizontal  level  flight  in  [13).  Using  this  approximation  and  a  backward 
difference  to  approximate  differentiation,  the  “discrete  time”  normal  acceleration  an(t) 
is 

/J4  J/(/)-2y(/-l)  +  j,(f-2)  _2 

a„(t)  = - y~2 - ms 

where  T,  seconds  is  the  sampling  interval.  In  order  to  penalize  excessive  normal  accel¬ 
eration,  Py  was  chosen  to  be  10(1  -  2:~'  +  c-2). 

The  same  cost  function  was  then  used  for  each  Example. 

5.1  Precompensator  Using  Current  Data 

Example  1.  Choosing  n,  =  0,  one  obtains  Sfc’1)  =  -0.1084.  The  dc  gain  of 
Y/Y'  is  then  0.850.  Simulated  performance  is  shown  in  Figure  2. 

Example  2.  With  n,  =  5, 

S(z~')  =  -0.0931  +  0.0093c’1  +  0.0033z“2  -  0.0074c’3  -  0.0039c  4  -  0.0173c’5. 

The  dc  gain  of  Y/Y‘  is  0.856,  Figure  3  shows  the  performance  obtained.  In  Figures  2 
and  3,  the  delay  between  terrain  peaks  and  aircraft  altitude  peaks  is  visible. 
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5.2 


Precompensator  Using  Future  Data 


Example  3.  With  n,  =  k  =  5, 

S(-'‘)  =  -0.0483  -  O.OOSOiT1  -  0.0090;'  2  -  0.007Gz  3  -  0.0019;-4  -  0.0363;"5. 

The  dc  gain  of  Y/Y ’  here  is  0.867.  Figure  1  shows  the  performance  obtained.  Here 
the  delay  between  terrain  peaks  and  aircraft  altitut.e  peaks  is  reduced  by  the  action  of 
the  precompensator  on  future  terrain  values. 

5.3  Preprogrammed  Control 

Example  4.  Figure  5  shows  the  simulation  performance  using  nv  =  630.  This  is 
close  to  the  best  possible  minimization  of  J.  The  plots  of  elevator  activity  and  normal 
acceleration  are  much  smoother  than  for  the  other  Figures. 

The  solution  to  eqn(15)  is  a  vector  of  631  values  of  v(t).  The  first  five  are: 

(v0  «i  n2  v3  v4)  =  (0.3173  0.5238  0.6357  0.7404  0.8400). 
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5.4  Comparison 

Table  1  contains  some  measures  of  system  performance  for  the  four  Examples. 
Table  1.  Comparative  performance  of  Examples. 


Ex 

599 

J  -  y(0l2 

<=0 

Umax 

U-min 

Uamax 

™a  min 

(y  y)max 

(y*  - 

1 

916103 

662j03 

5.01 

-4.73 

5.38 

-5.80 

110.5 

-79.5 

2 

870103 

7011O3 

3.46 

-3.72 

4.32 

-4.00 

111.8 

-92.6 

3 

540,03 

435103 

2.30 

-2.75 

2.99 

-2.50 

61.0 

-57.0 

4 

826102 

468102 

1.84 

-1.27 

1.43 

-2.02 

40.2 

-25.4 

ym  —  y  is  in  metres 

u  is  in  degrees 

n„  is  in  units  of  gravitational  acceleration  (g) 

Quantity  J  is  explicitly  minimized  in  these  Examples.  Note  that  the  use  of  pre¬ 
programmed  tracking  (Example  4)  gives  an  order  of  magnitude  improvement  in  perfor¬ 
mance  measured  by  the  quadratic  criteria  compared  with  the  precompensators.  This  is 
essentially  because  631  performance  enhancing  parameters  are  being  used,  rather  than 
one  or  six  as  in  the  precompensators. 

Quantity  E”oly*(<)  -  y(<)P  is  a  measure  of  tracking  accuracy  and  is  one  compo¬ 
nent  of  J  (neglecting  the  portion  after  t  =  600). 

Quantities  umaz  and  umm  are  the  extreme  values  of  elevator  deflection  during  the 
simulation  and  are  a  measure  of  control  activity.  Examples  3  and  4,  which  both  make 
use  of  future  terrain  data,  show  that  it  is  possible  to  obtain  improved  tracking  accuracy 
with  reduced  control  activity  compared  with  Examples  1  and  2  which  use  current  and 
past  terrain  values. 

The  peak  values,  namaz  and  naml„,  of  normal  acceleration  are  reduced  in  a  similar 
manner  to  the  peak  elevator  deflections  when  future  terrain  data  are  used. 
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The  quantity  (y*  —  y  )mnI  is  the  maximum  distance  the  aircraft  reaches  below  the 
desired  trajectory,  and  subtracts  from  available  ground  clearance.  Examples  3  and  4, 
which  use  future  reference  values,  give  significant  reduction  of  this  quantity  compared 
with  Examples  1  and  2,  which  do  not  use  future  data.  The  relatively  small  additional 
reduction  in  ground  clearance  obtained  using  preprogrammed  tracking  does  not  reflect 
the  order  of  magnitude  reduction  in  5^190|y*lO  —  y(t)]2  between  Examples  4  and  3. 

The  effect  of  varying  Pu(z~1)  and  Py(z~1)  has  not  been  shown  here.  For  a  given 
structure  of  v(t),  tracking  accuracy  may  be  improved  at  the  expense  of  increased  control 
activity  and  normal  accelerations  by  reducing  Pu  and  Py. 

6  FUTURE  DIRECTIONS 

For  implementation  in  a  more  practical  terrain  following  system,  the  algorithms 
would  need  to  be  modified  or  at  least  embedded  in  a  different  framework. 

The  requirement  of  prior  knowledge  of  the  entire  future  path  is  restrictive  with 
respect  to  both  mission  and  computational  requirements.  Several  proposed  methods 
for  avoiding  this  requirement  are  outlined. 

6.1  Deterministic  Reference  -  Precompensator 

A  simple  approach  is  to  design  the  precompensator  using  a  relatively  short  piece 
of  terrain  and  then  to  use  this  precompensator  over  the  whole  flight. 

An  extension  of  this  approach  would  be  to  partition  the  reference  signal  into  frames 
and  to  design  a  precompensator  for  each  frame.  Simulation  studies  would  be  necessary 
to  determine  appropriate  frame  lengths.  The  precompensator  design  assumes  zero 
initial  conditions  and  it  would  be  necessary  to  examine  the  effect  of  violating  this 
assumption  at  the  frame  boundaries. 

The  possibility  of  a  finite  horizon  cost  function  could  be  considered.  This  would 
mean  that  the  X  matrix  (with  dimension  n,  +  1)  would  not  have  an  autocorrelation 
structure  so  that  a  conventional  linear  equation  solver  would  be  required  for  the  solution 
of  eqn(10).  This  may  be  acceptable  if  the  precompensator  order  n,  is  fairly  small. 
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6.2  Deterministic  Reference  -  Preprogrammed  Tracking 

The  “patching”  together  of  frames  of  reference  trajectory  would  be  necessary  if 
preprogrammed  control  inputs  were  used  since  the  dimension  of  the  least  squares  prob¬ 
lem  (eqn(15))  is  approximately  equal  to  the  frame  length.  Studies  would  be  necessary 
to  determine  the  best  trade  off  between  the  good  performance  expected  with  a  long 
frame  length,  and  the  computational  burden.  Again,  zero  initial  conditions  are  assumed 
in  the  design  formulation  and  the  effects  of  violating  these  may  need  to  be  studied.  It 
is  possible  for  the  problem  to  be  ill-conditioned  and  this  may  need  to  be  examined. 

6.3  Stochastic  Reference  -  Precompensator 

The  reference  signal  could  be  represented  as  a  stochastic  process  obtained  from  the 
output  of  a  linear  dynamic  system  driven  by  white  noise.  Minimization  of  an  infinite 
horizon  cost  function  is  appropriate  for  this  formulation  because  of  the  relationship 
between  the  variance  of  the  output  of  a  system  driven  by  white  noise,  and  the  system’s 
impulse  response. 

System  identification  techniques  (probably  recursive  least  squares)  could  be  used 
to  identify  and  update  in  real  time  a  dynamic  model  of  the  terrain  using  a  frame 
of  future  terrain  values.  Studies  to  determine  the  appropriate  model  order  and  data 
discarding  strategies  would  need  to  be  carried  out. 

This  terrain  model  could  then  be  used  to  redesign  the  precompensator  as  the 
terrain  characteristics  varied.  Again,  computational  considerations  may  dictate  how 
often  this  could  occur.  The  precompensator  would  operate  on  actual  terrain  values. 

7  CONCLUSIONS 

Some  systematic,  optimization  based  approaches  for  improving  control  system 
tracking  performance  were  developed.  In  particular,  methods  for  incorporating  prior 
knowledge  of  the  reference  input  values  were  considered. 

The  described  approaches  were  applied  to  a  simplified  aircraft  terrain  following 
system  simulation.  It  was  demonstrated  that  use  of  prior  knowledge  of  terrain  height 
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values  could  give  improved  performance  in  the  sense  of  achieving  more  accurate  tracking 
with  less  control  activity.  This  allows  the  aircraft  to  fly  at  lower  altitude. 

Some  proposals  for  applying  the  algorithms  to  more  practical  terrain  following 
systems  were  presented. 
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FIGURE  1.  SCHEMATIC  DIAGRAM  OF  CON"ROt  SYSTEM 
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FIGURE  2.  PERFORMANCE  OBTAINED  USING  ns  =  0  AND  k  =  0  (EXAMPLE  1 ) 
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FIGURE  5.  PERFORMANCE  OBTAINED  USING  PREPROGRAMMED  TRACKING 
(EXAMPLE  4) 
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